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Abstract 

A precise calculation of the top quark pair production associated with a hard photon is essential 
for testing the electroweak property of the top quark in the Standard Model (SM). We investigate 
the one-loop QCD corrections to the process e~^e~ — t- tt'j at the International Linear Collider 
(ILC), and find that the iT-factor can be as large as 1.238 (1.105, 1.060) for a center-of-mass 
energy ^/s = 500 (800, 1500) GeV. The transverse momentum distributions of the top quark and 
photon are respectively shown at leading order (LO) and next-to-leading order (NLO). Due to the 
asymmetric rapidity distribution of the top (anti-top) quark, we also study the top quark forward- 
backward asymmetry (A^p^) in tt'j production at NLO, which is found to be 45.82 (55.25, 55. 
for = 500 (800, 1500) GeV. 

PACS numbers: 14.65.Ha,13.66.Bc,14.70.Bh 
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I. INTRODUCTION 



The Standard Model has been well testified by various experiments except that the 
Higgs scalar is still left as a missing piece. Therefore, testing the mechanism of electroweak 
symmetry breaking (EWSB) ^ is now an urgent task for particle physics. The top quark, 
which was discovered at the Tevatron in 1995 [sl, is distinguished for its large mass and 
short lifetime j4|. Since it is free from the QCD confinement, the top quark productions and 
decays are much cleaner than the light quarks ^j. Thus, the top quark is speculated to be a 
sensitive probe for the EWSB and new physics 6| . So far, except for the forward-backward 
asymmetry in ti production at the Tevatron [3], most of the measurements on the top quark 
are consistent with the Standard Model (SM) predictions, such as the cross sections of the 
ti production and the single top production |8|. 

Very recently, the CDF Collaboration has reported its observation of ^7 events with a 
luminosity of 6.0 fb~^ However, the small statistics still limits the precision study of 
the gauge coupling of ti'y. In addition, since the initial photon radiation severely affects 
the sensitivity of tt'y production to anomalous top quark couplings [l^, the SM electric 
charge of the top quark has not been measured directly up to now. Although the dominant 
contribution to the ti'-f production comes from the gluon fusion at the LHC, it is challenging 
to determine the top charge by measuring the cross section ratio a [ti'-f)/ a (tt) ll| due to the 
huge QCD backgrounds. The LO and NLO QCD calculations of the process pp{p) — > tt'y 
have been recently carried out at hadron colliders in Ref . . 

In contrast, as a clean top quark factory, the ILC will allow for a precison test for the 



top quark property 



13( 1 . Since the gauge coupling of tt jis sensitive to new physics, it has 



been studied intensively in a model independent way 



14l |. Some new physics models can 



jittle Higgs model 15|. It is also found that 



16| a sizable top quark electric (weak) dipole 



also affect the tt'y coupling sizably, such as the 
in the supersymmetric and multi-Higgs models 
momentum can be induced by the non-standard CP violating interactions. At the ILC, due 
to a high luminosity, such anomalous couplings of the top quark can be measured at the one 
percent level, which is much better than that at the LHC 13|]. Therefore, the high order 
calculations for the top quark processes at the ILC are needed to meet the experimental 
precision. The QCD and electroweak corrections to the process e"'"e~,77 — )■ ti have been 



studied in Ref. 



17j . The CP violation effects induced at loop level in the top pair production 



2 



through (un) polarized e'^e annihilation 18| or photon fusion 19| have been investigated 



in some extensions of the SM. In Ref. 



20|, the authors have studied the QCD corrections to 



the bb'j production at LEPl. In this paper, we calculate the one-loop QCD corrections to 
the process e~^e~ — )■ tt'-f at the ILC. 

This work is organized as follows. In section II, a brief description for the NLO QCD 
calculations is given. The discussions and numerical studies are presented in section III. 
Finally, the conclusions are drawn in section V. 

II. A DESCRIPTION OF ANALYTICAL CALCULATIONS 



In our calculations, the NLO QCD corrections {Aaqco) are divided into two parts: the 
virtual corrections (Aamr) and the real gluon radiation corrections (Aareai)- We adopt the 
dimensional regularization to isolate all the ultraviolet divergences (UV) in the one-loop 
amplitudes generated by FeynArts 21[ and remove them with the on-mass-shell renormal- 



ization scheme 



22] . The FormCalc-6.1 23| and LoopTools-2.5 2J] are employed to simplify 



the amplitudes and to perform the numerical calculations respectively. 

On the other hand, the infrared (IR) divergences arising from the contributions of virtual 
gluon exchange in loops are still left. According to the Kinoshita-Lee-Nauenberg (KLN) 



theorem 



25| . these IR divergences will be canceled by the real gluon bremsstrahlung cor- 



rections in the soft gluon limit. We denote the momentums of initial and final states for the 
real gluon emission process as follows: 



e+(pi) + e-(p2) -> t{ki) + t{k2) + 7(^:3) + 9{k) 



(1) 



We take the phase-space-slicing method to isolate the IR singularity and divide the real 
corrections into the hard and soft parts by the energy of the emitted gluon in the calculations 



Aareai = Aa soft + Aa hard ■ (2) 

we can obtain the soft part of the cross sections by the 



[261,1271: 

In the soft gluon approximation 
following equation: 

dAasoft = da^- 



28 



d'^k ( ki 



ko 



asCp 

27r2 Je,<ae, 2Eg \ki-k k2 ■ k 



(3) 
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where Cp = |, Eg = y\k\'^ + rn"^ and we give a fictitious mass nig to the gluon to ehminate 
the IR divergence. Note that this dependence on the non-physical mass will be exactly 
canceled by the virtual corrections which are also evaluated with a non-zero gluon mass. 
AEg is the energy cutoff of the soft gluon and we require Eg < AEg <^ The hard 

gluon [Eg > AEg) radiation corrections, which are insensitive to the small gluon mass, can 
be directly evaluated by the numerical Monte Carlo method 



Finally, the finite total cross section of the process e+e — tt'j including the LO (ao) and 
NLO QCD corrections {AaQco) can be expressed as 

o'tot = 0-0 + AaqcD = (70 + Aamr + Aasoft + Aahard = cro(l + Sqcd) , (4) 
where Sqcd = AaQcn/o'o is the relative QCD corrections at the order of 0{as). 

III. NUMERICAL RESULTS AND DISCUSSIONS 
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FIG. 1: (a) The dependence of the NLO QCD corrections on the soft cutoff bg in the process 
e+e~ — )■ tt7 for mg = 10~^ GeV at ^/s = SOOGeV; (b) The amplified curve marked with the 
calculation errors for AaqcD versus 6s- 



In the numerical evaluations, we take the input parameters of the SM as 30 1 



mt = 17L2 GeV, = 0.519991 MeV, mz = 91.19 GeV, 



sin^^w = 0.2228, a(m 



2 N-l 



127.918 



(5) 



Since the total cross section is independent of the non-physical parameters nig and the soft 
cutoff 5s{5s = AEg/ Ef,, Ef, = y/s/2), we display the curves of the NLO QCD corrections 
versus the cutoff 6s for rrig = 10~® GeV in Fig. 1(a), where we fix the renormalization scale 
/i = /io = mt. For the strong coupling constant as(/i), we use the two-loop evolution of 
it with the QCD parameter A"-'^^ = 226 MeV and get as(/io) = 0.1078. It can be seen 
that the values of Aatard and Aa^r + Acsoft vary with the change of the soft cutoff 6s, 
but the total NLO QCD correction Aagcn is not dependent on 6s within the reasonable 
calculation errors. In order to demonstrate this more clearly, we amplify the curve of AaqcD 
in Fig. 1(b). We also verify that the total correction is indeed independent on nig for the 
fixed 6s- Therefore, in the following calculations, we take 6s = 2x 10~^ and nig = 10~* GeV. 
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FIG. 2: The transverse momentum distributions of the top quark and photon at LO and NLO 
QCD respectively for the process e~^e~ — )■ ttj. The bands correspond to the variation of the 
renormahzation scale in the interval I-Iq/2 < jjL < 2fiQ. 

In Fig.2(a-b), we show the transverse momentum distributions of the top quark and 
photon at LO and NLO QCD for ^/s = 500 GeV. In the calculations, we maintain the 
electron mass and impose a transverse momentum cut > 15 GeV to exclude soft photon 
emission. It can be seen that the QCD corrections greatly enhance the magnitudes of the LO 
differential cross section dao/dpT- But the shapes of these distributions are not dramatically 
changed. Most of the top quarks are produced in the region of 40 GeV < < 150 GeV; 
while the photons are inclined to distribute in the region 15 GeV < < 60 GeV. We 
note that the detection of energetic photons produced by hard scattering goes through the 
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29.37(3) 


35.96(3) 


1.224(1) 


15 


21.62(2) 


26.76(2) 


1.238(1) 


20 


16.76(2) 


20.94(2) 


1.249(2) 


25 


13.37(1) 


16.87(1) 


1.262(1) 


30 


10.84(1) 


13.81(1) 


1.274(1) 



TABLE I: The LO cross sections, the NLO QCD total cross sections and i^-factors under different 
cuts for process e~^e~ — )■ ttj at t/s = 500GeV. 



definition of an isolation criterion [3l|. However, this relies on the detailed Monte Carlo 
simulation, such as parton shower, which is beyond the scope of our study. In table I, we 
present the effects of different p'^ cuts on the LO cross section, the NLO QCD total cross 
sections and the i^-factor of the ti'j production at a/s = 500 GeV. For a higher p'^ cut, the 
cross sections become smaller and the relative corrections get larger. 




FIG. 3: (a) The cross sections of e~^e~ — t- tt^ versus ^/s at LO and NLO QCD respectively; (b) The 
corresponding i^-factor versus ^/s. The bands correspond to the variation of the renormalization 
scale in the interval fio/2 < n < 2fj.Q. 

In Fig.3(a-b), we give the dependence of cross sections and relative corrections on the 
center-of-mass energy y/s. Since the process e~^e~ — )■ ti'j in the SM is induced by the pure 
electro- weak interaction at the order 0{a^), the LO cross section will not be affected by 
the variation of the renormalization scale in the strong coupling. However, the NLO QCD 
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27.30(2) 


1.263(1) 45.72(4) 
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55.3(2) 


14.41(3) 
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26.76(2) 


1.238(1) 45.82(4) 


27.06(3) 


1.105(1) 


55.4(2) 


14.34(3) 


1.060(2) 56.0(3) 


2/io 


26.31(2) 


1.217(1) 45.90(5) 


26.83(3) 


1.096(1) 


55.5(6) 


14.26(3) 


1.054(2) 56.1(3) 



TABLE II: The numerical results of the NLO QCD total cross sections, the iT-factor and the top 
quark forward-backward asymmetry at different values of the renormahzaiotn scale for the process 
e+e~ — )■ ti^. 

corrections are leading order in Us and show a weak dependence on the scale, due to the sup- 
pression of the loop factor. We display the values of the NLO QCD total cross section {(Ttot), 
the ii'-factor and the top quark forward-backward asymmetry at the renormalization scale 
/i = /io/2, /io, 2/io in table II. The uncertainty of the NLO scale dependence is approximately 
3.7% (1.8%,1.0%) for = 500 GeV (800 GeV, 1500 GeV) when the scale /i is varied between 
/io/2 and 2/io. The uncertainty is defined as 5 = [|cr(/io/2) — cr(/io)| + |(T(2/io) — cr(/to)|]/a"(/to). 
When setting the scale at /io, we find that the largest production rates of ti^ will reach about 
25.82 fb and 29.98 fb at LO and NLO QCD respectively around = 600 GeV, where the 
threshold effect may dominate. The corresponding relative QCD correction can be 16.1%. 
When a/s is greater than 600 GeV, the cross sections drop rapidly, due to the s-channel 
suppression. 

We also investigate the rapidity-differential cross sections of the top (anti-top) quark in 
Fig.4(a-b) and find that the events of the top (anti-top) quark for y > are more (less) 
than that for ?/ < 0. This asymmetry is caused by the huge interference effect between 



the photon and boson mediated in the process e^e~ — )■ ^* jZ^ — )■ tt'y 32[. The QCD 
corrections enhance the LO distributions of the top (anti-top) quark but do not distort 
their shapes significantly. In order to present this asymmetry, we can define the top quark 
forward-backward asymmetry (A^^) in the process e^e~ tt'y as 



N{y^ > 0) - N{y, < 0) 



^^^"iV(y,>0) + iV(y,<0)- 

Here N{yt > 0) and N{yt < 0) denote the events of top quarks moving along or against a 
given direction, which is chosen as the direction of the incoming particle e~ in our calcula- 
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FIG. 4: The rapidity distributions of top and anti-top at LO and NLO QCD respectively in the 
process e^e~ — )■ tt'^. The bands correspond to the variation of the renormahzation scale in the 
interval < < 2/io- 



tions. In Fig.5, we can see that the QCD corrections give arise to a negative contribution 
to the LO forward-backward asymmetry. We list the values of A^p^ at different scales in 
table. II. We find that the dependence of A*p^ on the renormahzation scale is very weak, due 
to the cancelation of strong coupling between numerator and denominator in Eq.(6). It is 
also noted that the value of this asymmetry is not sensitive to the collision energy when -y/i 
is greater than 900 GeV. The maximal values of A^p^ can reach 58.1% and 56.4% for /i = /ig 
at LO and NLO QCD respectively. 



IV. CONCLUSIONS 



In this paper, we discussed in detail the one-loop QCD corrections to the process 
e^e~ — )■ tt7 at the ILC. We found that the QCD corrections can significantly enhance 
the production rate of ti'-j and show a weak dependence on the renormahzation scale. The 
shapes of differential distributions of the top quark and photon are not be greatly affected 
by the QCD corrections. When fixing yU = /io, we found that the total cross section and the 
top quark forward-backward asymmetry can respectively reach 26.76 (27.06, 14.34) fh and 
45.82 (55.4, 56.0)% at NLO QCD for = 500 (800, 1500) GeV. 
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FIG. 5: The dependence of the forward-backward asymmetry Ap^ on -y/i in the production of 
tt^ at LO and NLO QCD respectively at the ILC. The band corresponds to the variation of the 
renormahzation scale in the interval hq/2 < fx < 2/io- 
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